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ABSTRACT 


Sound  waves  of  an  unknown  origin  have  been  observed 
to  radiate  from  specific  areas  in  the  Rocky  Mountains.  This  prompted 
early  researchers  to  refer  to  them  as  Mountain- Associated-Waves  or 
MAW.  These  waves  travel  at  sub-sonic  frequencies  (oscillation  periods 
1.0  to  1000  seconds)  and  are  detected  through  a  network  of  receivers 
located  thousands  of  kilometers  apart.  The  MAW  events  were  found 
to  be  associated  with  areas  of  maximum  wind  in  the  upper  atmosphere. 
This  suggests  a  production  mechanism  which  may  depend  on  the  speed 
of  the  air  flow  past  the  mountains. 

A  simple  dipole  model  has  been  used  to  investigate  the 
generation  of  MM.  Based  on  the  theory  of  airflow  around  cylinders, 
but  modified  by  incorporating  laboratory  results,  the  model  has  been 
extended  to  serve  on  a  synoptic  scale.  In  this  "natural"  dipole 
model  the  cylinder  has  been  replaced  by  the  massif  of  a  mountain. 

The  observed  MAW  events  were  found  to  have  preferred 
azimuthal  ranges  at  a  given  receiving  station.  This  may  be  due  in 
part  to  the  "steering"  influence  of  the  atmosphere  upon  sound  wave 
propagation. 


When  the  predicted  events  were  compared  to  observations, 
it  was  found  that  the  observed  events  were  too  small  in  number  to 
warrant  any  definite  conclusions.  However,  the  natural  dipole  appears 
to  be  a  possible  explanation  in  over  50  percent  of  the  cases  studied. 
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CHAPTER  1 


INTRODUCTION 


1.1  Purpose  of  Study 

The  low-frequency  band  of  of  the  acoustical  spectrum 
below  15  Hz  is  generally  referred  to  as  the  infrasonic  region. 

There  are  a  number  of  atmospheric  wave  structures  travelling  at 
acoustic  velocities  in  the  10  to  100  second  period  range.  Infra¬ 
sound  in  this  region  has  been  indentified  as  caning  from  man-made 
sources  such  as  nuclear  explosions  (Donn  and  Ewing,  1962)  and 
sonic  boons  (Goerke,  1971) .  Various  natural  sources  of  infrasound 
include  exploding  meteors  (Pekeris,  1948) ,  volcanic  explosions 
(Goerke,  et  al,  1965) ,  earthquakes  (Cook  and  Young,  1962) ,  certain 
auroral  events  (Campbell  and  Young,  1963;  Maeda  and  Young,  1966), 
and  severe  storms  (Bcwnan  and  Bedard,  1972). 

In  addition  to  these  infrasonic  emissions,  naturally 
occurring  infrasound  of  unknown  origin  is  frequently  detected. 

Work  by  Greene  and  Howard  (1975)  indicates  that  these  signals 
occur  mainly  in  winter  in  the  Northern  Hemisphere.  Location  of  a 
source  requires  triangulation  from  two  or  more  observatories.  The 
azimuth  range  for  each  station  receiving  these  unknown  signals 
shows  scare  predominant  bands.  There  is  a  clear  tendency  for  the 
source  locations  to  lie  in  rugged  mountainous  areas.  This  prompted 
early  investigators  to  refer  to  these  signals  as  "Mountain-Associated 
Waves"  (MAW) . 
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Fig.  1.1  shews  the  observed  distribution  of  azimuth 
angles  at  Edson,  Alberta.  The  predominant  signal  occurs  along 
azimuths  between  280°  and  320°,  suggesting  the  Rocky  Mountains 
for  possible  source  locations. 
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Fig.  1.1  Directions  from  which  Infrasonic 

Waves  are  received  at  Edson,  Alberta 

The  length  of  time  during  which  MAW  signals  are  con¬ 
tinuously  received  (often  called  an  "event")  at  a  station  varies 
from  hours  to  days.  Variations  in  signal  characteristics  during 
an  event  indicate  that  atmospheric  parameters  change  over  the 
course  of  observation,  the  nature  of  the  source  changes,  or  that 
two  or  more  sources  may  be  contributing  to  the  signal. 
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Craine,  Thomas  (1970A,  1970B,  1971) ,  and  co-workers  at 
Washington  State  University  and  the  University  of  Idaho  have  for 
several  years  investigated  natural  aerodynamic  sound  as  a  likely 
source  of  Mountain-Associated  Waves.  Larsen  (1971)  has  also  suggested 
the  interaction  of  winds  with  mountainous  regions  to  be  a  possible 
source  mechanism. 

1 . 2  Proposed  Theory 


This  thesis  will  attempt  to  show  that  natural  aero¬ 
dynamic  sound  of  dipole  origin  is  a  likely  source  for  many  of  the 
observed  MAW  events.  A  model  for  the  production  of  sound  at  infra- 
sonic  frequencies  in  the  atmosphere  will  be  constructed.  Propagation 
effects  will  be  discussed  to  point  out  the  difficulties  involved 
in  determining  the  origin  of  a  sound  ray  emitted  by  a  source  of 
uncertain  location. 

The  characteristics  (ie.  azimuth  of  arrival  and  intensity) 
of  the  sound  received  at  an  observatory  on  a  given  day  are  predicted. 

A  triangulation  region  is  constructed  based  on  the  model  and  then 
compared  to  actual  observations  for  the  same  time  period.  Conclusions 
are  then  drawn  according  to  the  correlation  between  predicted  and 


actual  measurements. 


CHAPTER  II 


AERODYNAMIC  SOUND 

2.1  Theory  of  Sound  Production 

Following  the  outline  of  acoustical  theory  and 
nomenclature  given  by  Bates  and  Stevens  (1966)  and  Kimpo  (1973) 
there  are  three  fundamental  ways  in  which  kinetic  energy  may 
be  converted  into  infrasonic  wave  motion,  namely  through  the 
action,  respectively,  of  monopole,  dipole  and  quadrupole  oscillators. 

2.1.1  Monopole 

A  monopole,  sometimes  called  also  a  simple  source,  is 
produced  by  forcing  the  mass  in  a  fixed  region  of  space  to  fluctuate. 

One  could  picture  a  pulsating  sphere  which  causes  density  fluctuations 
in  the  surrounding  fluid,  as  shown  schematically  in  Fig.  2.1a.  The 
density  fluctuations  then  propagate  radially  outward  in  the  form  of 
sound  waves,  producing  a  spherically  symmetrical  intensity  field. 

Since  the  inward  and  outward  fluctuations  involve  the  same  mass  flux, 
there  is  no  net  momentum  change. 

2.1.2  Dipole 

Forcing  the  momentum  in  a  fixed  region  of  space  to  fluctuate 
or,  equivalently,  forcing  the  rates  of  mass  flux  across  fixed  surfaces 
to  vary  constitutes  a  dipole  (See  Fig.  2.1b).  Fluctuation  of  momen¬ 
tum  corresponds  to  a  fluctuation  of  force  causing  a  density  variation 
in  the  surrounding  fluid  which  propagates  away  frcm  the  dipole  in  the 
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form  of  a  sound  wave.  Physically,  a  dipole  may  be  thought  of  as  two 
monopoles,  joined  by  a  fixed  axis,  pulsating  180  degrees  out  of  phase. 
The  intensity  of  the  sound  field  produced  by  a  dipole  is  a  maximum  along 
the  axis,  decreasing  to  zero  perpendicular  to  the  axis  (Fig.  2.1b). 

2.1.3  Quadrupole 


An  acoustic  quadrupole  is  equivalent  to  two  dipoles 
acting  together.  The  two  dipoles  can  be  thought  of  as  two  forces  a 
small  distance  apart  which  cause  a  fluctuating  shear  stress  within  a 
fluid.  The  intensity  field  produced  by  a  quadrupole  (see  Fig.  2.1c) 
has  maxima  along  lines  x  and  y,  and  minima  along  lines  xx  and  yy. 

The  monopole,  dipole  and  quadrupole  form  a  mono  tonic 
sequence  of  decreasing  efficiency  in  the  sense  that: 

(a)  a  dipole  is  a  less  efficient  means  of  energy 
conversion  than  a  monopole,  and  a  quadrupole  is  less 
efficient  than  a  dipole. 

(b)  the  differences  in  decreasing  efficiency  become 
more  pronounced  as  the  frequency  of  oscillation 


decreases . 


- 
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Fig.  2.1a 


X 


Fig.  2.1b 


Fig.  2.1c 


Fig.  2.1  Intensity  field  of  sound  produced  by  a)  monopole,  b)  dipole, 
and  c)  quadrupole . 
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2.2  Aerodynamic  Sound  Production 
2.2.1  Sound  Field  Characteristics 


The  production  of  sound  in  a  free  fluid,  ie.  a  fluid  not 
constrained  by,  nor  containing  within  it,  any  solid  boundaries,  will 
be  considered  next.  Sound  produced  in  the  absence  of  any  solid  boundaries 
is  called  aerodynamic  sound.  The  first  theory  which  concerned  itself 
with  the  intensity  of  sound  produced  by  a  fluctuating  fluid  flow  was 
developed  by  Ligh thill  (1951) . 

The  stresses  produced  in  the  fluid  give  rise  to  a  quadrupole 
type  radiation  field.  Although  Lighthill's  theory  assumes  that  there 
are  no  solid  boundaries  present,  he  pointed  out  that  they  might  well 
play  an  important  role  in  sound  production,  and  could  make  their 
presence  felt  in  two  ways: 

1)  Sound  generated  by  quadrupole s  will  be  reflected  and 
diffracted  by  solid  boundaries. 

2)  The  quadrupoles  will  be  no  longer  distributed  over  the 
whole  of  space,  but  only  throughout  the  region  external 
to  che  solid  boundaries. 

There  may  also  be  a  resultant  distribution  of  dipoles 
at  the  solid-fluid  interface.  (Dipoles  are  likely  since, 
in  acoustics,  they  correspond  co  excernally-applied 
fluctuating  forces,  and  such  forces  are  present  between 
the  fluid  and  solid  boundary) . 
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Curie  (1955)  extended  Lighthill's  theory  by  examining  the 
influence  of  solid  boundaries  upon  the  production  of  aerodynamic  sound. 

He  shoved  that  the  sound  field  may  be  regarded  as  being  derived  from 
twD  distinct  origins: 

1)  the  quadrupole  field  which  represents  the  fluctuating 
applied  stresses  within  the  fluid; 

2)  the  dipole  field,  which  represents  the  fluctuating 
force  with  which  the  solid  boundaries  act  on  the  fluid. 


The  relative  magnitude  of  dipole  and  quadripole  radiation 
intensities  is  given  as: 


-—a  — )  x  function  of  the  Reynolds  number 

Id  vc/ 


(2.1) 


Thus  for  a  sufficiently  small  Mach  number  (U/c) ,  the  contribution  to 
the  sound  field  frcm  the  dipoles  should  be  greater  than  from  the 
quadrupoles.  The  exact  value  of  Mach  number  for  this  to  occur  depends 
upon  the  nature  of  the  function  of  the  Reynolds  number,  which  in  turn 
is  determined  by  the  characteristics  of  the  flow. 


2.2.2  Flow  of  Air  Past  a  Cylinder 

A  stationary  cylinder  embedded  in  an  air  stream,  with  its 
axis  perpendicular  to  the  flow  may  be  expected  to  cause  quadrupole  and 
dipole  acoustical  fields  to  be  set  up.  Before  examining  the  characteristics 
of  the  sound  field  produced  in  such  a  manner,  it  is  useful  to  review 
the  concept  of  the  Reynolds  number. 


. 
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The  nature  of  a  fLow,  ie.  whether  it  is  laminar  or  turbulent, 
and  its  position  on  a  scale  indicating  the  relative  importance  of 
laminar  to  turbulent  tendencies  is  indicated  by  the  Reynolds  number, 
a  non-dimensional  quantity  defined  as: 


Re  = 


inertial  forces 
viscous  forces 


Ud p 

M 


(2.2) 


Defining  the  kinematic  viscosity  by  v 
may  be  written  as: 


M  ,  the  Reynolds  number 

P 


(2.3) 


2. 2. 2.1  Karman  vortices  and  Strouhal  Number 


The  first  extensive  experiments  on  the  flow  of  air  past  a 
stationary  circular  cylinder  were  done  by  Strouhal  (1878)  .  He  found 
that  sound  of  an  essentially  pure  tonal  character  is  radiated  for 
Reynolds  numbers  in  the  range  of  50  to  3  x  10^ .  Such  sounds  are 
usually  referred  to  as  "Aeolian  tones",  and  are  familiar  as  the 
"hunming"  of  telephone  or  power  lines  in  high  winds.  Strouhal  found 
that  the  frequency  of  the  note  was  predicted  by: 

f  =  0.185  V  (2.4) 

d 

It  is  now  customary  to  express  the  frequency  in  terms  of 
the  non-dimensional  Strouhal  number,  S,  such  that: 

S  =  LA  (2.5) 

U 

(S  =  0.185  in  Strouhal' s  experiments) 
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Since  Strouhal's  early  findings,  there  have  been  numerous  experiments 
dealing  with  flow  past  cylinders  in  an  attempt  to  explain  the  production 
of  aeolian  tones.  A  significant  advancement  occurred  when  Von  Karman 
(1911)  developed  his  now- famous  theory  of  the  "vortex  street".  The 
basic  vortex  theorem  states  that  the  fluid  flow  across  a  cylinder  will 
result  in  the  formation  of  two  unsymnetrical  rows  of  vortices  with 
alternating  arrangement  (Fig.  2.2) . 


Fig.  2.2  Production  of  Karman  vortices  in  the 
wake  of  a  cylinder 


According  to  this  theory,  an  arrangement  of  twD  rows  of 
vortices  is  unstable  unless  they  are  positioned  in  accordance  with  rhe 
spacing  ratio  given  by: 


Y  =  0.28 

A 


(2.6) 


This  phenomenon  of  vortex  shedding  is  responsible  for 
the  production  of  aeolian  tones;  the  frequency  of  the  tone  emitted 
is  equal  to  the  rate  at  which  the  vortices  are  shed. 
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Since  Von  Raman' s  original  work,  many  theoretical  and 
experimental  investigations  on  the  Karman  vorcex  street  have  been 
conducted  (eg.  Roshko,  1954a,  1954b  and  1961,  Gerrard,  1966) .  Although 
very  little  is  understood  about  the  principles  involved  in  the  actual 
vortex  shedding  at  the  cylinder,  some  interesting  results  have  been 
obtained.  For  example,  it  was  originally  thought  that  organized 
vortex  shedding  occurs  only  over  a  rather  limited  range  of  Reynolds 
numbers  (  50  <  Re  <10^  ).  Laboratory  results  have  shown,  hcwever,  that 
it  exists  to  Reynolds  numbers  of  at  least  10 ^  (Roshko,  1961)  and  there 
appears  to  be  no  reason  to  believe  chat  this  is  an  upper  limit.  Hcwever , 
it  is  extremely  difficult  to  produce  higher  Reynolds  numbers  in  a 
laboratory,  so  results  in  this  higher  range  remain  unavailable.  Another 
result  has  been  the  determination  of  the  frequency  with  which  the  vortices 
are  shed  from  the  cylinder,  and  hence  the  frequency  of  the  sound  produced. 

As  mentioned  earlier,  this  frequency  is  usually  expressed  in  terms  of 
a  Strouhal  number  (Eqn.  2.5).  The  Strouhal  shedding  frequency  appears  to 
be  about  0.28  for  circular  cylinders  at  high  Reynolds  numbers  (Roshko, 

1961),  decreasing  to  about  0.18  at  low  Reynolds  numbers  (Gerrard,  1955). 

It  should  be  noted  that  the  frequency  of  the  sound  produced 
by  the  shedding  of  vortices  is  independent  of  temperature,  density 
and  composition  of  the  fluid. 

For  circular  cones  the  frequency  is  determined  by  the  local 
diameter  and  there  is  a  spanwise  coupling  of  frequencies  along  the 
length  of  the  cone  (Gaster,  1969) .  Instead  of  observing  a  single 


' 
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frequency,  or  pure  tone,  as  in  the  case  for  cylinders,  one  could 
expect  a  broadening  of  the  spectrum.  This  effect  becomes  increasingly 
important  at  lower  frequencies. 

2.2.2. 2  Sound  Field 


Theory  shows  that  the  fluctuating  forces  of  an  air  flow  on 
a  rigid  body  may  be  expected  to  possess  an  associated  dipole  radiation 
field  (Lighthill,  1952) .  Experiments  have  confirmed  this  (Gerrard, 

1955)  .  The  sound  field  produced  by  an  air  flow  past  a  cylinder, 
perpendicular  to  the  axis,  is  found  to  possess  a  dipole  field  (Fig.  2.3) 
of  nearly  pure  tone  character;  the  harmonics  present  are  considerably 
less  intense  than  the  fundamental  frequency  (Etkin,  1957)  . 


Fig.  2.3  Resulting  dipole  field  formed  by  a  cylinder 
embedded  in  a  flow  of  air. 

The  sound  generated  is  aerodynamic,  ie.  it  is  not  produced 
by  vibration  of  the  cylinder  itself.  The  sound  spectrum  has  a  strong  peak 
at  the  fundamental  frequency,  a  weak  peak  ac  double  the  fundamental 
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frequency  and  wide-band  noise  at  higher  frequencies.  The  wide-band 
noise  appears  to  originate  from  quadrupole  sources  produced  in  the 
wake  of  the  cylinder.  The  fundamental  tone  radiates  most  strongly  in 
a  direction  perpendicular  to  the  flow,  while  the  weak  harmonic  radiates 
most  strongly  parallel  to  the  flow.  The  sound  produced  shews  no 
dependence  upon  the  physical  properties  of  the  cylinder  (eg.  type  of 
material,  hollow  vs.  solid  cylinder) . 

2.2. 2. 3  Intensity  Field 

Experiments  suggest  (Gerrard,  1955)  that  at  large  distances 
from  the  cylinder  axis  (£.>1)  in  an  air  flow  of  velocity  U,  relative 
to  the  cylinder,  the  intensity  of  the  sound  obeys  the  relation: 

I  =  U4  1 2 .9  cos2  Q  '  f  (R)  (2.7) 

c  yp- 


Fig.  2.4  Intensity  at  a  distance  r  from  the  cylinder 
axis 
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The  intensity  is  independent  of  the  elastic  properties  of 
the  cylinder  (ie.  cylinders  of  steel  and  brass,  as  well  as  hollow 
cylinders  produced  the  same  intensity  for  all  other  parameters 
remaining  constant) .  Holle  (1938)  investigated  the  dependence  of 
sound  intensity  upon  the  velocity  for  various  values  of  cylinder  diameter. 
He  found  no  dependence  on  diameter  except  at  high  velocities 
(  U  >  30  m/sec) . 


■ 


CHAPTER  III 


PRODUCTION  OF  MOUNTAIN-ASSOCIATED  WAVES 


3 . 1  Natural  Dipole  Generators 


The  flow  of  air  past  a  stationary  cylinder  may  be 
intuitively  extended,  on  a  global  scale,  to  the  flow  of  air  past  a 
prominent  mountain  peak,  with  an  accompanying  conversion  of  kinetic 
to  acoustic  energy.  A  dipole  radiation  field  will  be  produced  in 
similar  manner  co  that  produced  by  a  cylinder.  The  mountain-associated 
waves  can  now  be  thought  of  as  a  form  of  aeolian  tones  in  die  infrasonic 
frequency  range. 

Although  there  is  little  theoretical  or  experimental 
evidence  to  support  this  theory  for  the  production  of  mountain-associated 
waves,  soma  authors  have  suggested  such  a  means  of  production,  and  the 
need  for  work  in  this  area  (Larson,  1971,  Gossard  and  Hooke,  1975)  . 

The  Von  Karman  vortex  street  has  been  observed  in  the  wake  of  mountain 
peaks  with  the  aid  of  satellite  photography  (NASA,  1967)  . 

A  splendid  exanple  of  the  natural  occurrence  of  vortex 
shedding  on  a  synoptic  scale,  with  Reynolds  numbers  far  greater  than 
can  be  achieved  in  a  laboratory  (Re~1010)  ,  is  given  by  Teunissen  (1977). 
By  recons true cing  the  flow  pattern,  lb unis sen  obtains  a  Reynolds  number 
of  5.2  x  109  and  a  Strouhal  number  of  0.15.  Ihe  latter  value  agrees 
well  with  Strouhal  numbers  obtained  in  laboratory  experiments.  Recall 
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that  Von  Karman  concluded  the  vortex  street  would  remain  suable  only 

if  the  ratio  h/X  were  equal  to  0.28.  Teunissen  calculates  a  value 
of  h/X  to  be  0.3. 


Photographs  of  vortex  streets  on  a  synoptic  scale  have  become 
fairly  familiar  in  recent  years  (eg.  Bayliss,  1976)  .  The  importance 
for  the  existence  of  vortex  streets  on  synoptic  scales  comes  to  the 
fore  when  one  recalls  that  aeolian  tones  are  produced  by  the  breaking 
away  of  vortices.  Just  as  acoustical  aeolian  tones  are  detected,  for 
example,  as  the  humming  of  power  lines  in  the  wind,  infrasonic  aeolian 
tones  may  exist  as  M^W. 

Teunissen  (1977)  has  calculated  a  vortex  shedding  frequency 
of  l(T4/sec  for  the  flow  of  air  past  a  mountain  that  is  1300  m  high. 
Recalling  the  theory  states  that  the  sound  frequency  is  equal  to  the 
vortex  shedding  frequency,  this  implies  a  sonic  frequency  of  10”4  Hz  , 
certainly  well  into  the  infrasonic  range. 

If  mountain-associated  waves  are  actually  produced  by  a 
natural  dipole  it  would  be  a  main  source  since,  after  a  simple  source, 
the  dipole  is  the  next  most  efficient  method  of  converting  kinetic 
energy  to  acoustical  energy. 

3.1.1  Dipole  Model 

This  model  deals  with  one  type  of  natural  dipole  source 
in  the  atmosphere,  the  flow  of  air  past  a  mountain  peak.  Twd  conditions 
are  necessary:  a  mountain  of  large  diameter  projecting  well  above  the 


. 


surrounding  lesser  peaks,  and  an  appropriate  flow  of  air.  Rather 
than  working  with  every  mountain  in  a  given  region,  only  the  most 
prominent  peak  is  considered  in  the  analysis. 

There  are  ten  mountain  regions  in  the  area  under  study 
(Fig.  3.2  and  Table  3.1) .  The  resolution  possible  from  topographical 
maps  enabled  a  boundary  to  be  drawn  for  each  of  the  regions.  These 
source  locations  need  only  be  checked  for  wind  speed  and  direction 
to  provide  input  for  the  dipole  model. 


3. 1.1.1  Frequency  of  Mountain- Associated  Waves 


For  atmospheric  flow  past  a  mountain  the  Reynolds  number 
is  typically  around  10^.  For  very  high  Reynolds  numbers  such  as 
these,  the  corresponding  Strouhal  number  appears  to  be  about  0.27 
(Roshko,  1961) .  A  check  can  be  made  on  the  reasonability  of  the 
frequency  prediction  by  using  Strouhal' s  Law.  The  actual  frequencies 
are  known  (.01  Hz<f<.l  Hz) ,  S  is  known,  U  can  be  obtained  from  upper- 
air  charts,  d  can  be  derived  and  then  checked,  using  topographical 
maps,  as  being  of  reasonable  value.  For  typical  values  of  U  and 
f,  d  has  a  value  of  about  500  meters;  there  is,  therefore,  no  need 
to  perform  this  check  in  each  instance.  (Note  that  d  is  now  the 
diameter  of  a  mountain  peak) . 
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3. 1.1. 2  Intensity  of  Mountain- Associated  Waves 


Equation  2.7  introduced  in  the  previous  chapter  will 
be  used  to  predict  the  intensity  of  the  mountain-associated  waves 
at  a  receiving  station. 


The  function,  f(R),  has  only  been  obtained  for  Reynolds 
numbers  up  to  10^  (Gerrard,  1955) .  In  this  range  it  has  the  form: 

f (R)  =  3.6  x  109  (Re)3  (3.1) 

Since  there  has  been  little  experimental  work  done  at  higher 
Reynolds  numbers,  the  above  functional  relation  for  f(R)  will  be 
extended  and  applied  to  atmospheric  flows. 


Substituting  for  Re 


,  equation  2.7  becomes: 


I  = 


(3.6  x  109) 


d3  U7  y 2p 

3  2 

c  v  ° 


sin2  6 


(3.2) 


Fig.  3.1  Intensity  at  a  distance  r  from  a  source 
of  diameter  d  in  a  flow  of  velocity  U 
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An  explanation  of  the  variables  and  the  procedure  for 
obtaining  their  values  follows. 

The  peak  diameter,  d,  can  be  obtained  from  Strouhal's 
law,  or  by  use  of  a  topographical  map.  A  value  of  500  m  has  been  found 
to  be  a  reasonable  estimate  for  most  mountain  peaks. 

The  determination  of  the  physically  significant  height 
of  the  peak,  y,  poses  a  problem;  should  y  be  measured  from  sea  level, 
from  the  base  of  the  mountain,  or  some  other  reference  level.  Since  it 
is  the  interaction  of  an  isolated  peak  with  the  wind  that  is  important, 
the  height,  y,  is  taken  as  that  portion  of  the  peak  which  is  above  the 
average  height  of  the  surrounding  terrain.  A  visual  inspection  of  a 
topographical  map  shows  the  average  height  above  sea  level  of  the  terrain 
to  be  about  2500  meters  in  Regions  1  to  9  and  3000  meters  in  Region  10. 
(Peak  height,  and  height  above  sea  level  are  given  in  Table  3.1). 

The  distance  from  source  to  receiver,  r,  is  assumed  to  follow 
the  great  circle  arc  between  the  two. 

The  angle  6  and  wind  speed  U,  as  shown  on  Fig.  3.1,  are 
obtained  from  upper-air  charts. 

The  values  of  air  density,  p  ,  sound  velocity,  c,  and 
kinematic  viscosity  of  air,  v  f  are  calculated  from  meteorological 
observations  (see  Chapter  VI). 


■ 


The  intensities  at  a  receiving  station  from  each  of  the 
sources  can  be  ranked,  and  a  dominant  dipole  generator  picked  out. 

In  sunmary,  the  dipole  model  predicts  generation  of 
infrasound  whenever  air  flows  past  a  mountain  peak.  The  sound  will  be 
received  at  a  station  if  the  wind  veccor  U  and  radius  vector  r 
(see  Fig.  3.1)  are  not  parallel,  ie.  6^0  degrees.  The  sources  will 
be  ranked  according  to  intensity. 
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TABLE  3.1 


Prominent  Mountains  in  Each  of  the  Ten  Regions 


LOCATION 

MOUNTAIN 

HEIGHT  ABOVE 
SEA  LEVEL  (m) 

PEAK  HEIGHT 
(m) 

LAT. 

DEGREES 

NORTH 

LONG. 

DEGREES 

WEST 

1. 

Mt.  Robson 

3954 

1454 

53 

119 

2. 

Mt.  Sylvia 

2946 

446 

58 

124 

3. 

Keele  Peak 

2972 

472 

63 

131 

4. 

Mt.  Isto 

2761 

261 

69 

144 

5. 

Mt.  McKinley 

6194 

3694 

63 

151 

6. 

Mt.  Logan 

6050 

3550 

60 

140 

7. 

Mt.  Wadding  ton 

4042 

1542 

51 

125 

8. 

Mt.  Rainier 

4392 

1892 

47 

121 

9. 

Mt.  Whitney 

4418 

1918 

36 

118 

• 

o 
!— 1 

Mt.  Elbert 

4399 

1399 

i 

39 

106 
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Fig.  3.2  Major  mountain  areas  and  highest  peak  in 
each  area 


CHAPTER  IV 


PROPAGATION  OF  SOUND  IN  THE  ATMOSPHERE 

The  characteristics  and  state  of  the  atmosphere  have  a 
significant  effect  on  the  infrasonic  signal  as  it  propagates  from  source 
to  receiver.  It  may  well  be  that  propagation  effects  have  more  influence 
in  determining  the  qualities  of  a  signal  at  a  receiving  station  than 
the  variations  at  the  sources.  Although  much  is  not  yet  understood 
concerning  propagation  of  infrasound,  seme  of  the  better-known  effects 
are  described  in  this  chapter. 

4 . 1  Sound  Pressure 


As  a  sound  wave  passes  an  observer  it  causes  fluctuations 
in  atmospheric  pressure.  For  infrasound  of  natural  origin  these 
fluctuations  are  usually  in  the  range  of  0.1  to  100  dyne/ errr 
(Cook,  1969).  Atmospheric  pressure  at  sea  level  is  about  10  dyne/cm  . 

Infrasound  is  detected  by  a  series  of  microphones  which 
convert  the  pressure  fluctuations  to  variation  of  an  electric  current. 
The  passage  of  a  sound  wave  also  produces  small  oscillations  of  the 
air  particles  and  a  subsequent  variation  in  temperature.  Hot-wire 
anemometers  have  been  used  to  measure  infrasound  by  means  of  these 
tempera fi  rp  fluctuations.  Since  microphones  must  be  located  out-of- 
doors,  they  also  record  variations  not  due  to  infrasonic  waves.  For 
example,  turbulent  eddies  produce  temperature  and  pressure  variations. 
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IVbst  infrasonic  microphones  are  designed  to  respond  to  pressure  rather 
than  temperature  variations/  since  this  has  been  found  to  be  the  more 
reliable  method. 

4.2  Velocity  of  Infrasound 


The  speed  of  sound  in  air  is: 

c* 2  =  fP  %  T  (4*1) 

Thus,  the  speed  is  proportional  to  the  square  root  of  the  absolute 
temperature .  (For  air  at  a  temperature  of  20°  C  the  speed  of  sound 
is  about  344  m/sec.).  This  formula  is  applicable  for  all  sound 
waves  from  the  low  infrasonic  frequences  (f  =  0.001  Hz  )  to  the  high 
ultrasonic  frequences  (f  >  20, 000  Hz). 

The  literature  on  atmospheric  sound  usually  distinguishes 
four  different  velocities  and  defines  them  as:  (1)  phase  velocity; 

(2)  group  velocity;  (3)  signal  velocity,  and;  (4)  horizontal  trace 
velocity. 


The  phase  velocity,  c,  is  the  speed  at  which  a  surface  of 
constant:  phase  travels  chrough  a  medium,  ie.  the  atmospheric  air.  This 
velocity  is  a  function  of  temperature.  For  horizontal  propagation, 
where  the  temperature  is  relatively  constant  within  a  layer,  che  phase 
velocity  is  also  nearly  constant.  For  vertical  propagation,  where 
the  variation  of  temperature  is  large,  the  phase  velocity  is  a 
function  of  the  temperature  field. 


.. 


$  j . 
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The  group  velocity  is  defined  as: 

ca  =  ,  (4.2) 

g  d  k 

where  oj  =  2  irf  is  the  angular  frequency  and  k  =  2  tt/\  is  che  wave- 
number.  It  is  the  velocity  at  which  che  sound  energy  propagates  through 
che  medium. 


The  signal  velocity  is  defined  as: 

Cs  =  °  '  (4.3) 

where  D  is  che  distance  (from  source  to  receiver)  and  t  is  the  time 
of  transit.  In  general,  the  phase,  group  and  signal  velocities  are 
different  and  depend  on  the  structure  of  the  atmosphere. 

The  horizontal  trace  velocity,  c^  ,  is  the  speed  that  is 
actually  measured  at  an  observing  station,  since  infrasonic  microphones 
are  usually  located  on  the  Earth's  surface  and  in  the  same  horizontal 
plane.  This  velocity  depends  on  the  angle  of  incidence,  @  ,  of 

the  sound  ray  such  that  ch  =  c  s^n  P  (Fig*  4.1) . 
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4.3  The  Attenuation  of  Sound  in  the  Atmosphere 


There  are  four  main  factors  in  the  atmosphere  which 
contribute  to  the  absorption  of  infrasound:  viscosity,  heat  con¬ 
duction,  water  vapor,  and  relaxation  of  thermal  energy. 


The  attenuation  of  infrasound  in  the  atmosphere  due  to 
viscosity  and  heat  conduction  is  considerably  less  than  for  audible 
sounds  because  of  the  low  frequency  of  oscillation.  The  absorption 
coefficient  is  about  1.6  x  10“4  /  r2 B  where  B  is  the  baromstric 

pressure  in  dyne/an2.  For  a  plane  wave  in  the  lower  atmosphere 
(ie.  at  altitudes  below  90  km)  at  a  period  of  r  =  10  seconds,  the 
absorption  is  less  that  2  x  10”9  ^  (Cook,  1962) .  Hence,  the 
loss  due  to  this  effect  is  insignificant,  even  for  propagation  over 
distances  of  thousands  of  kilometers. 


Below  altitudes  of  10  km  the  absorption  due  to  water 

vapor  may  be  expected  to  be  significant.  The  exact  variation  of 

this  absorption  with  atmospheric  pressure  is  not  accurately  known  for 

infrasonic  frequencies.  It  is  estimated  to  be  as  large  as 

5  x  10”9  /r2  —  at  sea  level  (Gook,  1962).  Thus,  for  infrasonic 

km 

levels  at  r  =  10  seconds,  the  absorption  due  to  water  vapor  is  still 

q  dfo 

insignificant,  being  only  5  x  10"° 


At  very  low  frequencies  there  is  an  absorption  due  to 


relaxation  of  the  thermal  energy  stored  in  vibrations  of  the  diatonic 


•' 
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molecules  in  air.  This  is  estimated  to  be  about  10~6  ^  for  a  wave 
°f  7  =  10  seconds,  in  the  lover  atmosphere  (Cook,  1962) .  Again, 
this  is  a  relatively  snail  loss  and  of  little  importance  in  infrasonic 
transmission. 

The  absorption  of  infrasound  in  the  acmosphere  due  to  all 
four  effects  is  thus  small  enough  so  that  propagation  can  occur  over 
large  distances  without  any  substantial  loss  in  energy. 

4.4  Temperature  Effects 

Since  the  refraction  of  sound  waves  in  the  presence  of 
atmospheric  temperature  gradients  is  a  well-known  phenomenon,  only  a 
brief  summary  will  be  given  here  (See  Fleagle  and  Businger,  1963)  . 

The  atmosphere  may  be  imagined  to  be  composed  of  a  large 
number  of  quasi-horizontal  layers,  each  differing  in  temperature  from 
adjacent  layers.  Sound  waves  experience  refraction  at  the  boundaries 
betwaen  layers  (Fig.  4.2) . 


Fig.  4.2  Refraction  of  a  sound  ray  due  to  vertical 
atmospheric  temperature  gradient. 
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Snell's  law  applies  at  each  interface, 

sin  i^  sin  sin  i2 

Cl  “  c2  -  C2 


such  that: 


,  etc. 


(4.4) 


where  cn  is  the  phase  velocity  of  sound  in  layer  n.  Using  the 
expression  for  velocity  of  sound  (Eqn.  4.1)  and  assuming  that  cp 
is  constant  along  the  ray  gives: 


sin  i  =  constant  x  /r  (4.5) 

Hence,  a  decrease  of  temperature  with  height  causes  the  sound  ray  to 
bend  upward.  Conversely,  an  inversion  would  cause  a  downward  bending 
of  the  ray. 


To  evaluate  the  constant,  T  and  i  must  be  known  at 

soma  point  along  the  ray.  For  a  ray  which  becomes  horizontal  at  some 

point  (  i  =  90°  )  it  is  seen  that: 

sin  i  1  .  /  T  „ 

- - -  =  — —  ,  or  tan  l  =/  -  (4.6) 

/ t”  /T7"  V  T'-T 

where  T*  is  the  temperature  at  the  point  where  the  ray  is  horizon  cal. 
If  the  temperature  is  assumed  to  be  a  linear  function  of  height 
such  that  T  =  T*  -  7z  ,  where  T  represents  the  absolute  temperature 
at  any  reference  height  and  7  is  the  lapse  rate,  then,  with  tan  i 
replaced  by  dx  /  dz  : 


This  equation  may  be  integrated  numerically  to  calculate  the  path  of 
a  ray.  It  has  been  used  in  the  original  detection  of  the  stratospheric 


temperature  inversion. 


About  the  time  of  World  War  I  it  became  apparent  that 
the  sound  of  gunfire  could  be  heard  within  a  radius  of  about  100  km 
from  the  source  and  often  beyond  200  km,  but  not  at  distances 
between  100  and  200  km.  As  shown  in  Fig.  4.3,  within  the  central 
circle  the  wave  was  audible,  but  at  greater  radii  upward  refraction 
made  the  wave  inaudible.  However ,  at  still  greater  radii  beyond 
the  zone  of  silence  waves  were  refracted  downward  to  ground  level  by 
a  warm  layer  of  air  in  the  high  atmosphere.  It  is  now  recognized  that 
the  temperature  at  50  km  above  the  surface  in  mid- latitudes  may  be  on 
occasion  as  high  as  325°  K  (52°  C) . 


Fig.  4.3  (a)  Propagation  of  sound  from  an  explosion 

(b)  Audible  and  Inaudible  zones 


Fran  the  previous  discussion  it  is  clear  that  the 


infrasonic  signals  received  at  a  station  may  not  have  travelled 
in  a  straight  line  frcm  the  source.  Indeed,  while  one  station  may 
receive  signals  from  a  given  source,  another  station  at  a  somewhat 
different  location  may  not  receive  them. 


Infrasonic  waves  produced  by  a  natural  dipole  are 


refracted  away  from  the  earth  under  conditions  of  normal  atmospheric 
lapse  rate.  For  any  lapse  rate  there  is  a  critical  ray  which  strikes 
the  earth  horizontally;  beyond  this  distance  the  only  sound  reaching 
the  ground  is  due  to  scattering  and  usually  very  weak  (Fig.  4.4). 
Equation  4.7  may  be  used  to  determine  the  maximum  range  of  the  critical 
ray  in  the  case  of  a  constant  lapse  rate  and  constant  temperature, 
with  the  result: 


(4.8) 


For  a  lapse  rate  of  6.5^  C/km,  temperature  of  273®  K  and  a  mountain 
peak  some  5  km  above  the  mean  terrain,  any  infrasound  produced 
by  the  peak  should  be  undetectable  at  distances  beyond  30  km.  This 
would  correspond  to  the  first  zone  of  silence  indicated  in  Fig.  4.2. 


Thus  it  may  be  concluded  that  even  though  a  mountain 


peak  is  producing  infrasonic  waves,  they  may  be  undetectable  by  a 
nearby  receiving  station. 


- 
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Fig.  4.4  Path  of  infrasonic  ray  produced  by  a 
natural  dipole  at  height  z 


4.5  Wind  Effects 


Wind  effects  play  an  important  role  in  determining 
the  path  of  sound  waves.  Near  the  surface  of  the  earth,  where  the 
wind  speed  is  usually  only  a  very  small  fraction  of  the  speed  of 
sound,  such  effects  are  rarely  substantial.  Hovever,  at  the 
tropopause,  in  the  vicinity  of  the  jet  stream  at  about  12  km, 
the  wind  speed  can  be  as  much  as  25%  of  the  sound  speed.  This 
causes  dramatic  changes  in  the  course  of  a  sound  ray. 

Table  4.1  gives  a  brief  summary  of  the  winds  for  the 
area  over  the  continental  United  States. 


TABLE  4.1 


East-West  Zonal  Winds  at  Several  Geometric  Altitudes 
Altitude  Mean  values  of  wind  speed  m/sec 


km 

(+)  ' 

16  Oct.  to 

31  March 

50 

+52 

55 

+58 

60 

+61 

eastward,  (-)  -  westward 


1  April  to 

1  June  to 

16  Aug.  to 

31  May 

15  Aug. 

15  Oct. 

+2 

-42 

0 

-5 

-47 

+2 

-8 

-51 

+5 

(Cook,  1969). 
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The  mean  north-south  (meridional)  wind  speeds  are  less 
than  10  m/sec  at  any  time  of  the  year,  and  have  an  average  value  of 
6  m/sec  northward. 

From  Table  4.1,  it  is  evident  that  at  50  km,  in  the 
winter  months,  the  ground  speed  of  an  eastward  propagating  sound 
wave  will  be  enhanced  by  52  m/sec,  while  a  westward  propagating 
wave  will  be  retarded  by  52  m/sec  (relative  to  the  ground) .  The 
50  km  thick  atmospheric  layer  between  the  stratopause  and  the  surface 
of  the  earth  serves  as  a  waveguide  for  eastward  propagation  of  sound 
waves  during  the  winter  months,  but  not  for  westward  propagation. 

In  summer  the  situation  is  reversed;  the  waveguide  is  effective 
only  for  westward  propagation. 

The  component  of  the  wind  perpendicular  to  a  sound 
wave  will  cause  the  wave  to  be  deflected.  This  shift  or  azimuthal 
deviation  is  dependent  on  three  factors  (Rockway,  1972)*:  the 
length  of  the  propagation  path,  the  magnitude  of  the  wind,  and  the 
latitude  of  the  final  position  of  the  ray.  Rockway  found  that  a 
wave  travelling  frcm  45°  N  latitude  to  55°  N  latitude,  and 
experiencing  a  50  m/sec  wind  wauld  have  a  10  degree  azimuthal 
deviation  frcm  a  no-wind  case.  Obviously  this  could  be  important 
when  trying  to  determine  the  origin  of  a  sound  wave. 


*  See  Kimpo  1973. 
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Wind  shear  also  has  a  significant  effect  upon  the 
propagation  direction  of  a  sound  wave.  The  normal  situation  in  the 
atmosphere  is  one  of  increasing  wind  speed  with  height.  Waves 
moving  against  the  wind  will  be  refracted  upward,  while  waves 
moving  in  the  direction  of  the  wind  will  be  refracted  downward. 

(Fig.  4.5). 


Fig.  4.5  Effect  of  vertical  wind  shear  upon  upwind 
and  downwind  propagation  of  sound  waves. 


The  usual  effect  of  wind  shear  is  thac  the  distance 
from  the  source  in  which  infrasound  can  be  detected  is  decreased, 
in  the  direction  of  the  wind.  A  vertical  wind  shear  of  6  m/sec /km 
is  equivalent  in  refractive  pcwer  to  vertical  temperature  gradient 
of  10°  C/km. 


4.6  Summary  of  Propagation  Effects 


The  preceding  sections  show  that  a  sound  wave  under¬ 
goes  numerous  changes  enroute  from  source  to  receiver.  Wind  and 
temperature  place  the  main  role  in  determining  the  characteristics 
of  the  transmitted  sound  rays,  while  the  attenuation  of  the  infrasonic 
wave  is  negligible.  Allowing  for  the  effects  of  wind  and  temperature 
one  might  expect  to  predict  the  azimuth  and  elevation  angle  at 
which  a  sound  wave  would  be  received  from  a  given  source.  However, 
a  fairly  detailed  knowledge  of  the  wind  and  temperature  at  each 
point  along  the  way  is  required.  Even  if  these  data  were  available 
the  precise  computations  of  ray  geometry  would  be  difficult.  But, 
as  a  general  rule-of-thumb,  ray  azimuths  should  be  shifted  slightly 
to  the  east  in  the  winter  because  of  the  effect  of  the  prevailing 
westerly  winds  in  mid-latitudes.  Since  the  magnitude  of  this  shift 
on  a  given  day  is  not  known,  no  corrections  have  been  made  in  this 
study  for  azimuthal  variations  due  to  wind  or  temperature  effects. 

As  a  final  thought,  it  should  be  noted  that  the  azimuths 
from  which  the  greatest  percentage  of  sound  waves  are  received 
may  be  strongly  affected  by  the  properties  of  the  atmospheric 
waveguide  which  allows  only  certain  preferred  directions  of  propaga¬ 
tion.  Therefore  a  preferred  azimuth  range  does  not  necessarily 
imply  that  the  greatest  number  of  sources  are  actually  located  within 

this  range. 
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CHAPTER  V 


RECEPTION  AND  DIRECTION  FINDING 

5.1  Measurement  and  Recording 

The  system  used  to  detect  infrasonic  signals  at  each 
of  the  stations  is  discussed  by  Cook  (1969)  .  A  brief  summary  will 
be  presented  in  this  section. 

An  array  of  at  least  four  microphones  at  each  receiving 
station  determines  four  characteristics  of  infrasonic  waves  passing 
through  the  area:  1)  the  amplitude  and  waveform  of  the  incident 
sound  pressure;  2)  the  horizontal  phase  velocity;  3)  the  distri¬ 
bution  of  sound  wave  energy  at  various  frequencies,  and;  4)  the 
direction  of  propagation  of  the  wave. 

The  microphones  are  located  at  ground  level,  approximately 
in  the  same  plane,  and  about  7  km  apart. 

Effects  of  pressure  fluctuations  due  to  local  turbulent 
eddies  are  minimized  by  a  noise  reducing  system  described  by  Daniels 
(1959) ,  which  consists  of  lines  of  pipe,  having  capillary  inlets, 
and  attached  to  the  microphone. 

The  received  signals  are  usually  recorded  on  paper 
charts,  one  at  each  of  the  four  microphone  sites.  The  characteristics 
of  the  sound  wave  are  then  obtained  by  overlaying  pairs  of  recordings 
from  the  four  charts.  This  method  requires  an  accurate  timing  trace 
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for  each  record  and  also  assumes  that  the  sound  waves  have  approxi¬ 
mately  plane  wavefronts. 

5 . 2  Direction  Finding 

By  overlaying  the  chares  and  visually  cor re lacing 
the  waveforms  from  pairs  of  microphones  che  azimuth  of  the  approach 
of  the  wave  and  the  horizontal  trace  velocity  can  be  obtained  from 
a  comparison  of  the  time  of  arrival  of  the  waveform  at  each  micro¬ 
phone.  Detection  of  signals  is  often  difficult  since  local  noise 
may  obscure  the  desired  signal,  and  superposition  of  signals  arriving 
from  two  or  more  directions  may  cause  the  waveform  to  vary  at  each 
microphone . 


Since  it  is  often  difficult  to  separate  mountain-associated 
wave  characteristics  from  outside  influences,  an  "event"  is  assumed 
eo  occur  when  a  narrow  range  of  azimuths  is  observed  for  an  ex¬ 
tended  period  of  time  (usually  greater  than  2  hours) .  A  variation 
in  azimuth  my  be  due  to  several  factors.  As  mentioned,  variation 
in  waveforms  ac  each  microphone  can  cause  poor  correlations.  Charts 
can  be  overlayed  with  an  accuracy  of  -  0.5  seconds  for  excellent 
signals.  The  quality  of  a  signal  is  determined  by  the  azimuth 
scatter  obtained  by  overlaying  pairs  of  charts.  Excellent  signals 
show  a  Jg  to  1  degree  scatter,  while  poor  signals  show  a  6  to  10 
degree  scatter.  Signals  of  greater  scatter  are  ignored  (Larson,  1971) . 

As  noted  in  the  previous  chapter  the  path  of  a  sound 
ray  may  be  altered  nurrerous  times  enroute  from  source  to  receiver. 


i  •  _ 


The  azimuth  angle,  in  particular,  varies  according  to  the  strength 
of  the  cross-wind  conpcnent.  Diamond  (1963)  has  shown  that  the 
azimuthal  error  due  to  a  cross-wind  velocity  can  be  estimated  as: 

ERROR  =  tan”1  J  (5.1) 

where  v  =  cross-wind  velocity 
c  =  speed  of  sound 

« 

Taking  the  possible  direction-finding  errors  into 
account,  Larson  (1971)  states  "...  azimuth  angles  for  good  signals 
do  not  have  errors  greater  than  10  degrees." 

5.2.1  Single  Source 


The  most  direct  path  of  a  sound  ray  is  along  the  great 
circle  route  between  a  source  and  a  receiver.  If  no  distortions 
of  the  ray's  path  were  present,  there  would  be  a  single  azimuth 
angle  of  arrival  at  a  station.  This  is  not  the  case,  howaver. 

Ins cead  a  range  of  azimuths  should  be  expected.  If  only  one  source 
is  predicted  by  the  dipole  model,  then  there  would  be  a  10  degree 
error  in  calculating  the  azimuth.  In  other  words,  a  range  of 
azimuths  of  -  5  degrees  around  the  true  heading  is  the  best  that  can 
be  expected  in  most  cases  (Fig.  5.1). 


■ 
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Fig.  5.1  Resolvable  azimuth  range  around  the 
true  heading. 


5.2.2  Two  or  More  Sources 


If  the  directions  of  two  or  more  sources  are  to  be 
determined  by  the  dipole  model,  then  an  azimuth  error  of  -  5  degrees 
must  be  expected  in  each  of  the  bearings  of  the  respective  sources. 
This  means  that  the  azimuths  of  the  true  headings  of  two  sources 
must  be  at  least  10  degrees  apart  for  resolution  between  the  two 
to  be  possible  (Fig.  5.2). 

In  practice,  however,  direction  finding  must  at  times 
be  somewhat  subjective.  For  example,  consider  the  situation  in  Fig. 


5.3. 


' 
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Fig.  5.2  Predicted  azimuth  ranges  at  a  receiving 
station  from  three  sources  a,  b,  and  c. 


1.  a  i>  10°  ~  resolution  between  sources  a  and  b 

(or  c)  is  possible 

o:  2  <  10°  -  resolution  between  sources  b  and  c 
is  not  possible 


2. 


' 
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Fig.  5.3  Azimuth  range  for  a  receiver  at  Edson  (ML) 
for  sources  at  mountain  locations  1,7  and  8 

In  this  case  the  dipole  model  has  predicted  that  a  receiver  at 
Edson,  Alberta  (ML)  should  receive  three  distinct  signals;  one 
signal  from  mountain  location  number  1,  a  second  signal  from  location 
7,  and  a  third  signal  frcm  location  8.  The  true  bearings  from 
sources  1  and  7  differ  by  only  8  degrees;  thus  it  should  not 
be  possible  to  resolve  signals  produced  by  these  two  sources.  On 
the  other  hand,  since  source  8  has  a  bearing  differing  by  some  44 
degrees  frcm  that  of  source  7,  one  could  assume  that  these  two 
sources  my  be  easily  distinguished.  However,  since  the  terrain 


between  the  three  sources  is  very  mountainous,  and  all  three  locations 
are  producing  signals,  one  may  suspect  that  signals  may  also  be 


emitted  by  mountains  in  the  region  between  sources  7  and 
resolution  between  any  of  these  three  mountain  locations 
to  be  very  difficult.  An  azimuth  range  as  shown  in  Fig. 


8.  Thus, 
appears 
5.3  is 


assumed. 
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CHAPTER  VI 


COMPARISON  BETWEEN  PREDICTED  AND 
OBSERVED  MOUNTAIN-ASSOCIATED  WAVE  CHARACTERISTICS 


6.1  Infrasonic  Triangula tion  Events 


The  first  step  in  this  analysis  of  infrasonic  signals 
from  unknown  sources  is  the  determination  of  probable  triangulation 
events.  A  triangulation  event  is  defined  as  an  interval  of  time 
during  which  infrasonic  signals  are  received  at  the  triangulating 
stations.  A  natural  requirement  is  than  the  stations  receive  signals 
at  the  appropriate  relative  times.  Events  are  selected  for  study  on 
the  basis  of  the  following  criteria: 

1)  The  events  cover  a  relatively  long  interval  of  time, 
preferably  greater  than  2  hours.  This  improves  the 
chances  that  a  single  dominant  source  is  being  observed. 

2)  The  triangulation  area  is  stable,  ie.  the  azimuth 
angles  are  fairly  steady  throughout  the  observation 
interval;  this  again  increases  the  chances  that  a  single 
event  or  stable  configuration  of  events  is  being  detected. 

3)  Meteorological  data  must  be  available  for  the  time 
during  which  the  event  was  recorded. 
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Four  stations  vere  involved  in  the  detection  of  infrasonic 


signals  and  the  subsequent  determination  of  triangulation  points: 


1.  Fairbanks  (or  College) ,  Alaska 

2.  Pul lman-Washington ,  Washington 

3.  Boulder,  Colorado 

4.  Edson,  Alberta 


64. 51°N,  147 - 43CW 
46. 44°N,  117. 10°W 
40. 01°N,  105. 17°W 
53. 35°N,  116. 26°W 


(Edson  was  a  temporary  station  oaken  out  of  service 
in  1973.) 


Following  the  format  adopted  by  ocher  authors,  the  following 
abbreviations  will  be  employed  for  convenience: 

Edson  -  ML  Pullman  -  PW 

Fairbanks  -  CO  Boulder  -  BO 

(See  Fig.  6.1) 

In  order  to  locate  a  triangulation  area,  a  signal  must  be  received  by 
at  least  two  of  these  stations  at  the  proper  relative  times. 

The  location  and  times  of  the  triangulacion  areas  which 
were  chosen  for  study  are  given  in  Table  6.1  (also  Fig.  6.1). 
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TABLE  6.1 

Infrasonic  Triangulation  Events  for  Aug.  1971  to  Feb.  1972 

EVENT  DATE  TIME  LATITUDE  LONGITUDE 

(Greenwich  Mean)  (Degrees  North)  (Degrees  Vfest) 


1 

30/09/71 

0030-0800 

55 

126 

2 

05/10/71 

0000-0100 

53 

125 

3 

07/10/71 

2130-2300 

46 

122 

4  a) 

08/10/71 

0000-1130 

56 

125 

b) 

08/10/71 

2300-2400 

60 

130 

c) 

09/10/71 

0000-0130 

59 

129 

5 

18/11/71 

0800-1600 

57 

128 

6 

01/12/71 

0400-0600 

60 

155 

7 

29/12/71 

0030-0200 

57 

130 

8 

14/01/72 

0030-0230 

46 

114 

9 

18/01/72 

0730-1330 

41 

no 

10 

25/01/72 

0200-2300 

51 

125 

11 

26/01/72 

0000-2400 

54 

122 

12 

27/01/72 

0000-0500 

54 

123 

13 

01/02/72 

0030  0930 

60 

120 

14 

02/02/72 

0030-0530 

57 

125 

15 

03/02/72 

0030-0430 

56 

120 

16 

09/02/72 

0730-1200 

58 

127 
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Fig.  6.1  Infrasonic  triangulation  events  for 
Aug.  1971  to  Feb.  1972. 
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6 . 2  Meteorological  Data 

The  dipole  model  requires,  as  input,  the  wind  speed  and 
direction  at  each  mountain  location.  It  is  also  desirable  that  these 
data  be  observed  near  the  peak  of  the  mountain. 

The  data  are  obtained  from  upper  level  charts  which 
contain  such  information  as  wind  speed,  wind  direction,  contour 
height,  temperature,  etc.  at  different  levels  in  the  atmosphere. 

Since  measurements  are  rarely  taken  directly  at  a  mountain 
location,  the  wind  velocity  at  the  ten  source  sites  must  be  estimated 
using  neighboring  data.  A  brief  explanation  of  the  method  used  to 
obtain  the  wind  velocity  at  each  location  follows. 

6.2.1  Wind  Velocity 


The  measurement  of  wind  velocity  is  needed  near  the 
peak  of  the  mountain.  The  available  upper-air  maps  were  the  700-mb 
charts  which  correspond  to  a  height  of  about  3000  m  above  sea  level. 
In  order  to  estimate  the  wind  speed,  the  geostrophic  wind  equation 
is  used: 


where 

and 


I  &A  (6.1) 

f  <5H 

f  =  2il  sin  <p  is  the  Goriolis  parameter 

v  =  geostrophic  wind 
y 

g  =  gravitational  acceleration  (9.81  m/sec^) 


« 


. 
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/I  =  Earth's  angular  velocity  (0.729  x  10  Vsec) 

<p  =  latitude  at  which  the  geos  trophic  wind  is 
to  be  calculated 

<5Z  =  standard  contour  interval  (60  gpm  or  approx. 

60  m) 

<5H  =  horizontal  spacing  of  contour  lines 

Using  a  latitude  of  0  =  50°  N  gives  a  value  for  tie 
geos  trophic  wind  as: 

Vg  =  5.23  x  106/  <$H  (6.2) 

from  which  Vg  merely  becomes  a  function  of  the  spacing  of  the  contour 
lines. 


The  primary  assumptions  underlying  Equation  6.1  are  that 
the  frictional  force  and  the  acceleration  are  small  in  comparison  to 
the  pressure  and  Coriolis  forces.  A  strict  balance  between  pressure 
and  Goriolis  forces  is  thus  implied  which  is  true  only  if  the  contour 
lines  are  straight  and  parallel.  Since  most  contours  are  curved,  an 

error  of  up  to  10%  may  be  introduced.  Also,  the  Coriolis  parameter 

varies  with  latitude.  However,  the  error  produced  by  assuming  no 
variation  is  negligible  in  mid-latitude .  The  geostrophic  equation 
also  ignores  terrain  effects  by  assuming  that  the  frictional  forces 
are  small.  At  low  altitudes  this  would  be  a  serious  source  of  error; 

however,  at  700-mb  the  frictional  forces  are  very  small  and  can 

usually  be  ignored.  The  geostrophic  wind  Equation  (6.1)  will,  in 
most  cases,  provide  an  overestimate  of  the  actual  wind  speed. 


'I 


. 
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The  wind  speed  at  the  various  mountain  peaks  can  be 
obtained  merely  by  measuring  the  spacing  betv^en  the  contour  lines 
and  applying  Equation  6.2.  The  wind  direction  is  also  easily 
obtainable  since  the  geostrophic  assumption  implies  that  it  be 
parallel  to  the  contour  lines. 

In  practice  it  may  not  be  necessary  co  use  Equation 

6.2  to  estimate  the  wind  speed.  By  taking  note  of  the  speed  ac 
neighboring  stations  a  good  estimate  can  often  be  made  at  the  source 
location. 

6 . 3  Computational  Procedure 


The  intensity  of  the  dipole  radiation  is  given  by 
Equation  3.1.  For  the  accuracy  attainable  in  the  calculation  of 
the  intensities  it  is  sufficient  to  assume  a  constant  temperature  of 
-15°C  (258°  K)  at  700-mb  over  the  period  in  which  the  infrasonic 
events  were  recorded  (August  to  February) .  The  air  density  and 
sound  velocity  (which  are  functions  of  temperature)  will  have  the 
values : 

p  =1.3  kg/m 3 
c  =320  m/sec 


—5  2 

The  kinematic  viscosicy  of  air  is  1.2  x  10  m  /sec  and 
the  average  peak  diameter  was  estimated  to  be  500  m.  Inserting  these 
values  into  Equation  3.1  gives: 


U7  y2  .  2q 
— sin  6 


I 
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(6.3) 


< 


Using  a  700-mb  upper- air  chart  corresponding  to  the 
period  of  observation  of  an  infrasonic  event,  the  following  procedure 
is  performed  for  each  of  the  ten  source  locations: 

1)  Estimate  the  wind  speed,  U. 

2)  Calculate  the  angle,  0  ,  for  each  station. 

3)  Obtain  the  peak  height,  y,  and  distance  to  station, 
r  (Table  6.2). 

4)  Use  the  above  values  for  U,  0  ,  y ,  and  r  to  predict 
the  dipole  intensity  at  each  of  the  four  stations. 

Once  the  intensities  frcm  each  source  at  a  given  station 
are  calculated  they  are  ranked  according  to  order  of  magnitude.  This 
serves  to  separate  the  dominant  source,  or  sources,  for  each  station. 
The  dominant  source  is  assumed  to  have  an  intensity  of  one  or  more 
orders  of  magnitude  greater  than  the  remaining  sources. 

Frcm  this  information  a  range  of  azimuths  at  each  station 
is  predicted  as  outlined  in  Chapter  V.  The  intersection  of  the 
azimuths  frcm  two  or  more  of  the  receiving  stations  defines  a  triangu¬ 
lation  region.  This  region  can  then  be  compared  to  the  observed 
triangulation  point  for  the  infrasonic  event  under  study. 


•• 


TABLE  6.2 


Peak  Heights  and  Distances  From  Receivers  to  Sources 


Mountain  Location 

Peak  Height,  y 

(m)  Distance  to  Receiver, 

r  (km) 

ML 

GO 

FW 

BO 

1 

1454 

186 

2049 

741 

1783 

2 

446 

729 

1411 

1378 

2427 

3 

472 

1410 

806 

2082 

3111 

4 

261 

2247 

521 

2918 

3945 

5 

3694 

2262 

243 

2792 

3933 

6 

3550 

1657 

569 

2176 

3318 

7 

1542 

648 

1957 

793 

1980 

8 

1892 

845 

2544 

334 

1499 

1918 

1399 


1895  3681  1125  1205 

1759  3874  1203  141 
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6.4  Results 


The  results  are  interpreted  mainly  on  the  degree  of 
proximity  between  the  predicted  region  and  the  observed  triangu¬ 
lation  point.  In  this  way  a  subjective  judgement  is  made  as  to 
the  likelihood  of  a  dipole  being  the  source  for  the  observed  MAW 
event. 


A. 


The  results  of  the  dipole  model  are  given  in  Appendix 


Since  I  «  IT7  in  the  dipole  model,  the  predicted 
triangulation  region  could  be  expected  to  be  in  the  areas  of 
maximum  wind  speed  ie.  the  jet  stream.  This  was  indeed  true  in  the 
majority  of  cases.  A  plot  of  the  latitude  of  the  observed  MAW 
event  versus  the  latitude  of  the  500-mb  jet  stream  is  shown  in 
Fig.  6.2.  The  correlation  between  the  two  suggests  that  the  wind 
speed  may  be  a  major  factor  in  the  production  of  MAW. 
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Fig.  6.2  Latitude  of  MAW  Event  and  Jet  Stream  for  each  of 
tiie  16  cases 


Examination  of  the  results  from  the  dipole  model  contained 
in  Appendix  A  indicates  that  when  closely  spaced  contour  gradients 
(ie.  strong  winds)  were  present  over  the  mountains,  both  the  pre¬ 
dicted  triangulation  area  and  the  observed  source  point  were  usually 
located  in  the  region  of  maximum  wind.  This  is  particularly  evident 
in  cases  2,  4,  5,  7,  8,  9,  15,  and  16.  Under  the  influence  of  a 
tight  gradient  the  predicted  areas  tended  to  follow  a  fairly  organized 
pattern.  When  the  gradient  was  slack,  however,  the  predicted  regions 
were  ill-defined  and  showed  very  little  correlation  to  the  observed 
triangulation  point.  This  is  evident  from  cases  1,  11,  12,  13, 
and  14. 


Cases  3  and  5  shew  a  well-developed  gradient  across 
the  mountains,  but  the  observed  and  predicted  regions  do  not  coincide. 
In  particular.  Case  3  shows  the  predicted  region  in  southern  Alaska, 
an  area  with  strong  upper  winds.  However,  the  observed  triangulation 
point  is  in  Washington,  in  the  center  of  a  high  pressure  area  with 
only  light  winds. 

Case  5  occurred  while  strong  winds  prevailed  over  much 
of  southern  Alaska  and  the  British  Columbia  coast.  The  predicted 
triangulation  regions  coincide  with  the  zone  of  maximum  winds.  The 
existence  of  three  separate  areas  is  most  likely  a  result  of  sub¬ 
dividing  the  mountainous  regions  into  three  distinct  source  locations. 
It  is  thus  likely  that  the  dipole  model  could  predict  a  region  ex¬ 
tending  along  the  coast  from  southern  British  Columbia  to  the  Gulf 
of  Alaska.  This  could  mean  that,  though  the  observed  point  falls 
within  the  boundary  of  this  predicted  region,  there  is  such  great 


. 
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uncertainty  in  the  precise  location  of  the  predicted  area  that  it 
would  be  difficult  to  deduce  any  useful  correlations. 

Case  6  also  shows  the  actual  triangulation  point  to  be 
well- removed  from  the  axis  of  maximum  wind.  The  dipole  model  does 
predict  a  source  at  the  observed  point;  however,  the  major  source 
is  located  in  the  tight  gradient  area  well  to  the  south. 

Case  4  is  particularly  interesting  in  that  it  involves 
three  distinct  triangulation  points  located  in  a  tight  contour 
gradient  area.  The  predicted  region  is  also  within  this  tight  gradient 
and  actually  encompasses  two  of  the  three  points,  the  third  being 
only  slightly  outside  this  area. 

Additional  testing  of  the  model  was  done  by  considering 
cumulative  statistics  for  the  16  cases.  The  predicted  arrival 
azimuths  from  which  signals  are  received  for  each  station  are  shown 
in  Fig.  6.3  through  6.6.  The  predominant  signal  emitters  for  the 
receiver  at  Edson  are  observed  at  an  azimuth  of  260  degrees.  There 
are  three  lesser  azimuths  of  308,  252  and  208  degrees  (Fig.  6.3). 

This  may  be  compared  to  Fig.  1.1  which  shows  the  majority  of  signals 
arriving  from  about  300  degrees.  The  College  data  predict  the  major 
azimuth  at  138  degrees  (Fig.  6.4).  The  Pullman-Moscow  observatory 
has  forecast  azimuths  of  273  degrees  and  324  degrees  (Fig.  6.5). 

Boulder  shows  a  preferred  azimuth  at  222  degrees  and  another  lesser 
heading  from  325  degrees  (Fig.  6.6). 


. 
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Since  observed  arrival  azimuth  statistics  are  available 


only  for  the  Edson  receiver,  no  comparison  was  made  between  predicted 
and  observed  results  for  the  other  three  receivers. 
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Fig.  6.3  Percentage  of  signals  received  along 

predicted  arrival  azimuths  at  Edson,  Alberta 


360 


Fig.  6.4  Percentage  of  signals  received  along 

predicted  arrival  azimuths  at  College,  Alaska 
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Fig.  6.5  Percentage  of  signals  received  along 

predicted  arrival  azimuths  at  Pul lman-Moscow, 
Washington 


360 


Fig.  6.6  Percentage  of  signals  received  along  predicted 
arrival  azimuths  at  Boulder,  Colorado 


It  is  customary  to  measure  the  power  of  a  signal  rather 
that  the  intensity.  According  to  Larsen  (1971)  the  observed 
pressure  amplitudes  for  MAW  vary  from  0.5  /zb  to  about  7.5  /zb 
zero- to- peak  (as  compared  to  typical  atmospheric  pressures  of  1000  mb) . 
However,  no  comparisons  have  been  made  between  observed  and  theoretical 
sound  pressures  in  the  present  study,  since  the  simple  model  does 
not  take  into  account  the  effects  of  attenuation  and  anomalous 
propagation. 


CHAPTER  VII 


SUMMARY  AND  CONCLUSIONS 


Theoretical  and  experimental  findings  have  shown 
that  the  flew  of  air  around  a  stationary  cylinder  will  result  in 
the  conversion  of  kinetic  to  acoustical  energy  for  sufficiently 
large  Reynolds  numbers.  The  sound  produced  in  this  manner  is  termed 
aerodynamic  sound  in  that  it  is  not  a  result  of  vibrations  of  the 
cylinder  itself,  but  rather  accompanies  the  shedding  of  vortices 
in  the  wake  of  the  cylinder.  The  associated  sound  field  possesses 
a  dipole  pattern  with  the  main  axis  perpendicular  to  the  flow. 

By  enlarging  the  concept  of  a  flew  of  air  past  a 
cylinder  to  an  analogous  flow  of  air  past  a  prominent  mountain  peak 
a  model  for  the  generation  of  Moun tain- Associated- Inf rasonic-Waves 
has  been  developed.  The  observation  of  Karman  vortex  streets  in 
the  lee  of  mountain  peaks  added  further  support  for  the  theory  of 
a  natural  dipole  generator. 

In  Chapter  V  the  propagation  of  infrasound  through  the 
atmosphere  was  examined.  It  is  evident  that  the  path  of  a  sound 
wave  may  vary  greatly  due  to  changes  in  the  state  of  the  atmosphere. 
Since  ray-tracing  theory  requires  a  detailed  knowledge  of  winds  and 
temperature,  it  was  considered  beyond  the  scope  of  this  study  to 
correct  propagation  paths  by  allowing  for  the  effects  of  varying 
wind  and  temperature  conditions.  The  main  conclusion  reached  is 
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that  the  largest  percentage  of  sources  may  not  actually  lie  within 
a  preferred  range  of  azimuths,  since  propagation  effects  may  determine 
the  direction  from  which  signals  can  be  detected. 

Moreover,  sensing  and  recording  of  infrasonic  signals 
includes  errors  which  make  direction  finding  difficult;  hence,  a 
range  of  azimuths  should  be  expected  rather  than  a  single  heading. 

In  most  cases  the  intersection  of  azimuthal  ranges  from  two  or 
more  of  the  receiving  stations  defined  a  triangulation  region. 

Wind  data  corresponding  to  the  times  of  observed  MAW 
events  were  used  in  the  model  to  predict  signal  source  regions. 

The  results  were  judged  according  to  the  proximity  of  the  observed 
and  predicted  areas. 

In  9  out  of  16  cases  examined  (case  numbers  2,  4,  5, 

6,  7,  8,  9,  15,  16)  it  appears  that  a  dipole  oscillator  could  be 
responsible  for  the  production  of  the  observed  MAW  events.  However, 
one  cannot  assume  that  it  is  the  only  mechanism  generating  the  ob¬ 
served  signal.  When  a  well-developed  contour  gradient  exists  across 
the  mountains ,  the  dipole  model  predicts  regions  which  are  in  fairly 
close  proximity  to  the  observed  points,  but  when  the  flew  is  weak 
the  predictive  skill  of  the  model  is  low. 

It  is  important  to  keep  in  mind  that  the  input  for  the 
model  cones  from  data  collected  at  one  specific  time  of  the  day, 
whereas  a  MAW  event  may  occur  over  a  period  of  several  hours. 
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Changing  wind  conditions  during  the  day  may  cause  a  corresponding 
change  in  forecast  triangulation  areas.  It  may  be  that  the  upper- 
air  data,  reported  at  0000  Z  and  1200  Z,  do  not  adequately  represent 
the  wind  pattern  prevailing  during  the  time  of  the  recorded  M£W 
event. 


The  700-mb  upper-air  chart  was  used  to  obtain  wind  data. 

This  chart  corresponds  to  an  altitude  of  approximately  3000  M  above 
sea  level.  The  highest  peaks,  such  as  Mt.  Robson,  Mt.  Fairweather 
and  Mt.  Logan,  range  frcm  about  4000  to  over  6000  m  above  sea  level, 
well  above  the  700-mb  level.  The  wind  conditions  will  normally  be 
different  at  the  level  of  these  peaks  because,  in  most  instances,  the 
wind  speed  increases  with  height.  If  calculations  for  the  higher  peaks 
are  based  only  on  700-mb  winds,  then  the  sound  intensity  would  be 

7 

expected  to  be  too  low.  Since  I  ct  U  ,  the  intensity  should  increase 
with  height.  The  higher  peaks  may  thus  be  still  dominant  dipole  sources 
and  the  resultant  triangulation  region  need  not  be  affected  to  any 
great  extent.  The  direction  of  the  wind  has  only  a  slight  effect  on  the 
intensity  of  infrasound  produced  by  the  dipole  model,  unless  the  wind 
vector  is  very  nearly  parallel  to  the  radius  vector  frcm  the  source  to 
the  receiver.  Thus,  even  if  the  direction  of  the  wind  changes  with 
height,  this  effect  will  most  likely  be  overshadowed  by  a  change  in  the 
speed  of  the  wind. 

Observation  shews  that  the  majority  of  the  signals 
received  at  Edson,  Alberta  arrive  from  about  300  degrees  (Fig.  1.1) . 

The  predicted  azimuths  of  arrival  for  Edson  (Fig.  6.3)  do  shew  that 
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a  large  percentage  of  signals  could  be  expected  to  ccme  fran  this 
heading;  however,  the  majority  cone  frcm  about  260  degrees.  The 
source  producing  these  signals  is  most  likely  Mt.  Robson,  which  is 
only  sane  180  km  from  the  Ed  son  receiving  array.  As  outlined  in 
Chapter  IV  on  propagation,  a  receiver  situated  100  to  200  km  frcm 
a  source  may  be  unable  to  detect  signals.  If  this  were  the  case 
with  Edson  and  Mt.  Robson,  then  a  number  of  signals  caning  from  the 
direction  of  260  degrees  may  be  undetectable,  meaning  the  largest 
percentage  of  signals  would  arrive  frcm  about  310  degrees. 

From  a  sample  of  only  16  cases  it  is  difficult  to 
draw  any  definite  conclusions  concerning  the  degree  to  which  the 
dipole  model  accounts  for  the  observed  MAW  events.  Results  frcm  the 
model  show  good  agreement  with  observations  in  a  number  of  cases, 
suggesting  that  dipole  oscillators  may  generate  Mountain- Associated- 
Waves  under  high  wind  conditions.  Other  mechanisms  for  the  production 
of  infrasound  are  likely  contributing  as  well.  "Natural"  quadrupoles, 
formed  in  the  turbulent  flow  over  the  mountains  as  well  as  in  clear- 
air- turbulence,  may  also  be  important  sources  of  infrasonic  signals. 

To  date,  little  experimental  wcrk  has  been  carried  out 

o 

at  very  high  Reynolds  numbers  (  >  10  )  corresponding  to  atmospheric 
flow  around  obstacles  the  size  of  mountains.  It  may  be  that  the 
function,  f(R),  used  in  calculating  intensities  is  different,  at 
such  high  Reynolds  numbers,  from  experimental  results  obtained  in 
the  laboratory. 


- 

g  j'ld  dy-  *  i.«r"  >  «  ■' 

■ 


Since  Mountain- Associated-Waves  are  a  fairly  recent 
discovery,  more  data  are  needed  before  any  definite  conclusions  can 
be  drawn  as  to  their  origin  and  production  mechanism. 


67 


REFERENCES 


Bates,  A.E.  and  R.B.  Stephens,  1966:  Acoustics  and  Vibrational  Physics. 
Edward  Arnold  Ltd.,  London.  Chapter  15. 

Bayliss,  P. ,  1976:  Weather,  31,  346. 

Bowman,  G.G. ,  S.  Howard  and  A.J.  Bedard,  1971:  Observations  of  Infra- 
sound  Subsonic  and  Disturbances  Related  to  Severe  Weather. 

Geophys.  J.R.  Astron.  Soc.  26,  215-242. 

Campbell,  W.H.  and  J.M.  Young,  1963:  Auroral- Zone  Observance  of  Infrasonic 
Waves  Related  to  Ionospheric  Disturbances  and  Geomagnetic  Activity. 

J.  Geophys.  Fes.,  6Q_,  5909-5916. 

Gook,  R.K. ,  1969:  Atmospheric  Sound  Propagation.  Atmospheric  Explora¬ 
tion  by  Remote  Probes,  National  Academy  of  Sciences  Panel  on 
Remote  Atmospheric  Probing,  2,  633-669. 

Gook,  R.K. ,  1962:  Strange  Sounds  in  the  Atmosphere.  Sound  1 (2) ,  12-16. 

Craine,  L.B.  and  J.E.  Thomas,  1970a:  Infrasonic  Waves  Observed  from 
Mountainous  Regions.  Paper  delivered  at  annual  meeting  of  Am. 

Geophys.  Union,  Washington,  D.C.  (April,  1970) . 

Craine,  L.B.  and  J.E.  Thomas,  1970b:  Infrasonic  Waves  in  the  Pacific 
Northwest.  Continued  proposal  resulting  in  National  Science 
Foundation  Grants  GA-23176  and  GA-23203  (June,  1970) . 

Curie,  N. ,  1955:  The  influence  of  Solid  Boundaries  upon  Aerodynamic 
Sound.  Proc.  Roy.  Soc.  London,  A  231,  505-514. 

Daniels,  F.B.,  1959:  Noise-Reducing  Line  Microphone  for  Frequencies 
Belcw  1  cps.  J.  Acoust.  Soc.  Am.,  31,  529-531. 

Diamond,  M. ,  1963:  Cross-Wind  Effect  on  Sound  Propagation.  J.  Appl. 
Meteorol.,  3,  208-210. 

Donn,  W.L.  and  M.  Ewing,  1962:  Atmospheric  Waves  from  Nuclear  Explosions- 
Part  II,  The  Soviet  Test  of  October  30,  1961.  J.  Atmos.  Sci. , 

19,  264-273. 

Etkin,  B. ,  G.K.  Korbacher  and  R.T.  Keefe,  1957:  Acoustic  Radiation 
from  a  Stationary  Cylinder  in  a  Fluid  Stream.  J.  Acoust.  Soc. 

Am.,  29,  30-36. 


Fehr,  U. ,  1967:  Measurements  of  Infrasound  from  Artificial  and 
Natural  Sources.  J.  Geophys.  Res. ,  72,  2403-2417. 

Fleagle,  R.G.  and  J.A.  Businger,  1963:  An  Introduction  to  Atmospheric 
Physics.  Academic  Press,  New  York.  Chapter  7. 

Gaster,  M. ,  1969:  Vortex  Shedding  from  Slender  Cones  at  Low  Reynolds 
Number.  J.  Fluid  Mech. ,  38,  565-576. 


68 


Gazaryan,  Y.L. ,  1961:  Infrasonic  Normal  Modes  in  the  Atmosphere. 

Soviet  Physics-Acoustics,  1_,  17-22. 

Gerrard,  J.H. ,  1963:  The  calculation  of  the  Fluctuating  Lift  on  a 
Circular  Cylinder  and  its  Application  to  the  Determination  of 
Aeolian  Tone  Intensity.  Advisory  Group  for  Aeronautical  Res. 
and  Dev.  (AGARD) ,  Report  463. 

Gerrard,  J.H. ,  1955:  Measurements  of  the  Sound  frcm  Circular  Cylinders 
in  an  Air  Stream.  Proc.  Phys.  Soc.  London,  68,  453-461. 

Goerke,  V.H.  and  M.W.  Woodward,  1966:  Infrasonic  Observation  of  a 
Severe  Weather  System.  Monthly  Weather  Review,  94,  395-398. 

Goerke,  V.H. ,  J.M.  Young  and  R.K.  COok,  1965:  Infrasonic  observations 
of  the  May  16,  1963  Volcanic  Explosion  on  the  Island  of  Bali. 

J.  Geophys.  Rev.,  70 (24) ,  6017-6022. 

Gossard,  E.E.  and  W.H.  Hooke,  1976:  Waves  in  the  Atmosphere.  Scientific 
Publishing  Co.,  New  York.  Chapter  9. 

Greene,  G.E.  and  J.  Howard,  1975:  Natural  Infrasound:  A  One  Year 
Global  Study. 

Holle,  Von  W. ,  1938:  Akust.  Z.,  3^,  321. 

Hunphreys,  W.J. ,  1964:  Physics  of  the  Air.  Dover  Publications, 

New  York,  Chapters  2  and  3. 

Kimpo,  F.L.,  1973:  Natural  Aerodynamic  Sound  as  a  Source  of  Mountain- 
Associated-  Infrasonic  Waves.  Master's  Thesis,  Dept,  of  Electrical 
Engineering,  Washington  State  Olivers ity. 

Larsen,  R.J. ,  1971:  Correlation  of  Wind  and  Geographical  Features 

with  the  Production  of  Certain  Infrasonic  Signals  in  the  Atmosphere. 
Geophys.  J.  Roy.  Astron.  Soc.,  26,  201-214. 

Lighthill,  M.J.,  1962:  On  Sound  Generated  Aerodynamical ly.  Proc. 

Roy.  Soc.  London,  A267,  147-182. 

Lighthill,  M.J. ,  1952:  On  Sound  Generated  Aerodynamically,  I,  General 
Theory.  Proc.  Roy.  Soc.  London,  A211,  564-587. 

Lighthill,  M.J.,  1954:  On  Sound  Generated  Aerodynamically,  II, 

Turbulence  as  a  Source  of  Sound.  Proc.  Roy.  Soc.  London,  A222, 

1-32. 

Maeda,  K.  and  J.  Young,  1966:  Propagation  of  Pressure  Waves  Produced 
by  Aurora.  J.  Geomagn.  Geoelectr. ,  Kyoto,  18,  275-299. 

IVfeecham,  W.C. ,  1971b:  On  Aerodynamic  Infrasound.  J.  Atmos.  Terres. 

Phys. ,  33,  149-155. 

NASA,  1967:  Earth  Photographs  from  Geminii  III,  IV,  and  V,  123. 


•' ' 


. 


Perkeris,  C.L.,  1948:  The  Propagation  of  a  Pressure  Pulse  in  the 
Atmosphere.  Phys.  Rev.,  73,  145-154. 

Phillips,  O.M. ,  1956:  The  Intensity  of  Aeolian  Tones.  J.  Fluid 
Mech.  1,  607-624. 

Rand  McNally,  1969:  The  International  Atlas,  Rand  McNally  and  Co. , 
New  York,  166-173. 

Roshko,  A.,  1953:  On  the  Development  of  Turbulent  Wakes  From  Vortex 
Streets.  National  Advisory  Conmittee  on  Aeronautics  (NACA) . 

Tech.  Note  2913,  1-77. 

Roshko,  A. ,  1954a:  On  the  Development  of  Turbulent  Wakes  from  Vortex 
Streets.  NACA  report  1191,  1-20. 

Roshko,  A. ,  1954b:  On  the  Drag  and  Shedding  Frequency  of  Two- 
Dimensional  Bluff  Bodies.  NACA  Tech.  Note  3169,  1-29. 

Roshko,  A.,  1961:  Experiments  on  the  Flow  Past  a  Circular  Cylinder 
at  Very  High  Reynolds  Number.  J.  Fluid  Mech.,  10,  345-356. 

Strouhal,  V.,  1878:  V.  Wied.  Ann.,  _5,  216. 

Teunissen,  H.W. ,  1977:  Global  Flow-Visualization.  Atmosphere,  15, 
54-56. 

Thomas,  J.E.  and  L.B.  Craine,  1971:  Acoustic  Gravity  Wave  Propagation 
in  a  Measured  Atmosphere.  Paper  submitted  for  publication  to 
Geophys.  J.  Roy.  Astron.  Soc. 

World  Meteorological  Organization,  1967:  The  Airflow  over  Mountains. 
WMO  Tech.  Note  #34,  Geneva,  Switz. ,  1-133. 


■ 


70 


APPENDIX  A 

PREDICTED  TRIANGULATION  REGIONS 

This  appendix  contains  the  results  obtained  frcm 
the  dipole  model.  The  shaded  areas  indicate  the  predicted 
triangulation  region  for  the  production  of  I^foun tain- Associated- 
Waves.  The  observed  triangulation  points  are  indicated  by  scars. 

Immediately  following  is  an  example  of  the  method 
used  to  predict  the  triangulation  region.  Case  2  is  selected 
since  it  is  reasonably  straight  forward.  Table  A1  contains  the 
intensities  expected  at  the  receivers  frcm  each  of  the  mountain 
locations  as  obtained  from  the  model,  by  use  of  Eqn.  6.3.  The 
Edson  (ML)  receiver  shows  two  dominant  sources  in  mountain 
locations  1  and  7  (ie.  Mt.  Robson  and  Mt.  Waddington  respectively). 
The  College  (CO)  receiver  has  dominant  sources  in  mountain  locations 
6  and  7  (Mt.  Logan  and  Mt.  Waddington  respectively) ;  Pullman  (PW) 
has  dominant  sources  in  locations  1  and  7;  Boulder  (BO)  has  a 
dominant  source  in  location  7  only.  Using  Edson  as  an  example,  a 
line  is  drawn  from  Source  1  to  the  receiver  at  Edson  (Fig.  Al) . 

Using  direction  finding  methods,  an  azimuth  range  is  assumed.  The 
sane  is  done  for  Source  7  to  Edson.  In  this  way  a  complete  azimuth 
range  is  constructed  for  infrasonic  signals  arriving  at  the  receiver 
at  Edson.  The  sane  procedure  is  followed  at  each  of  the  remaining 
three  stations.  Finally  a  triangulation  region  is  defined  by  the 
intersection  of  azimuth  ranges  from  each  receiver  (shown  as  the 


' 


. 

■ 

. 
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shaded  area  in  Fig.  Al)  . 


For  the  sake  of  clarity  only  the  predicted  triangulation 
region  is  shown  for  the  remaining  cases. 


TABLE  A1 


Predicted  Intensities  at  receiving  stations 
(in  units  of  nt  nT^  sec-^  x  10”^) 


Mountain  Station 

Location 


ML 

CO 

PW 

BO 

1 

*6160 

29 

*529 

30 

2 

2 

1 

4 

0 

3 

0 

0 

0 

0 

4 

0 

0 

0 

0 

5 

0 

3 

0 

0 

6 

51 

*115 

15 

9 

7 

*1030 

*941 

*3810 

*430 

8 

0 

0 

2 

0 

9 

2 

0 

5 

0 

10 

0 

0 

0 

0 

*  -  most  likely  source  fron  wliich  signals  would  be  received 
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Cose  2  700  mb  OOZ  05  Oct  71 


Fig.  A1  Method  used  to  construct  the 
predicted  triangulation  area. 
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Case  2 


700  mb  001  05  Oct  71 


Case  3 


700  mb  OOZ  08  Oct  71 
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Case  4 


700  mb  OOZ  09  Oct  71 


78 


700  mb  OOZ  18  Nov  71 


Case  5 


Case  6 


700  mb  OOZ  01  Dec  71 


Case  7 


700  mb  OOZ  29  Dec  71 
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Case  8 


700  mb  OOZ  14  Jan  72 


82 


120 «  110°  100 ° 


70° 


Case  9 


700  mb  UZ  18  Jan  72 


Case  10 


700  mb  121  25  Jan  72 
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Case  11  700  mb  OOZ  26  Jan  72 


85 


Case  12  700  mb  001  27  Jan  72 
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Case  13  700  mb  OOZ  01  Feb  72 


. 


2940 


87 


120°  no °  ioo° 


Case  14 


700  mb  OOZ  02  Feb  72 
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Case  15  700  mb  OOZ  03  Feb  72 


Case  16  700  mb  12Z  09  Feb  72 


